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The paper contains generalization of the renormgroup algorithm for boundary value problems 



in 
o 

. of mathematical physics and related concept of the renormgroup symmetry, formulated earlier by 

i authors with reference to models based on differential equations. These algorithm and symmetry 

are formulated now for models with non local (integral) equations. We discuss in detail and 
illustrate by examples applications of the generalized algorithm to models with not local terms 



which appear as linear functionals of the solution. 

1 Introduction 



< 
(N 

•n I The concept of RenormGroup Symmetry (RGS) appeared in mathematical physics in the 

^ I beginning of 90s [I1I21 (see, also reviews ^E]) being borrowed from theoretical physics. 

OO ■ In turn, the concept of the Renormalization Group (RG) for the first time has appeared 

^ ■ in the most complex branch of last — in the Quantum Field Theory (QFT). A presence 

O . of a group structure (Lie groups of transformations) in the QFT calculation results was 

^ ! discovered in the beginning of 50s by Stueckelberg and Petermann [5] (see, also [EIIZI)- 

^1 This structure and exact symmetry of a solution underlying have been then used by 

' N.N.Bogoliubov for creation of regular method for improving of an approximate solution 

^ ■ of QFT problems - the Renormalization Group Method (=RGM) IB] (see also [HIIIIIIIII]). 

This method contains elements of the theory of continuous group transformations (theory 
of Sophus Lie). Improvement of approximation properties appears to be most essential in 
^ ; a vicinity of the solution singularity. 

' Turn now to the fact that in the case of problems described by complicated equations 

^ — as, e.g., in the transfer theory (integro-differential Boltzmann equation) or in the 

quantum field theory (an infinite chain of engaged integro-differential Dyson- Schwinger 
equations) — only some components of solution or their integrated characteristics obey 
enough simple symmetry. So, in QFT the central object in RG transformations is the so 
called "function of an invariant charge" a (or "running coupling constant"), representing 
specific product of Lorentz-invariant amplitudes of propagators di , vertices and the 
expansion parameter a 

a = aT^ Y\ di . 

i 

At the same time, Schwinger- Dyson equations include only functions di Tk and the cou- 
pling constant a , separately, but not their product d . In the one-velocity plane transfer 
problem the RG-invariance property is related to the asymptotics of "density of particles. 
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moving deep into the medium" n^{x) , a; — *■ cxo not entering the Bohzmann equation^. 

The renormahzation group concept was transferred to mathematical physics ^112] with 
the same pragmatic goal, as in QFT, in mind — "improvement" of the solution behaviour 
in vicinity of a singularity. Remind that proliferation of the RG method from QFT to 
others fields of theoretical physics (the theory of critical phenomena, physics of polymers 
and so on) has caused various and sometimes essentially different (in comparison with ini- 
tial) forms of realization of RG ideology (see, for example, the review |12jJ). In application 
to problems of mathematical physics there appeared different variants of formulations of 
RG method [H-EIl- 

For the boundary value problems (BVPs) of mathematical physics using DE, we devel- 
oped (see, e.g., [1112], and also recent reviews W, p. 232-249], [22j the RG algorithm, which 
is quite distinct from earlier ones. To this difference became obvious, we remind, that in 
a basis of the RG Method, as it was initially formulated by N.N.Bogoliubov and one of 
authors jH] for QFT problems, lays the use of an exact group property of a solution. One 
of the well-known formulations of this property is the functional equation (representing 
only a group composition law) for the invariant charge in QFT. In every concrete case, 
revealing of similar symmetry (i.e., of group property) for the solution demands a special, 
usually non-standard, analysis (see, for example, discussion in the papers [2S1 1121 CSl), 
that is an algorithmic drawback of the RG technique. 

Coming back to mathematical physics (MP) we note, that here we usually deal with 
the problems based on systems of DEs the symmetry of which can be found by a regular 
way with the help of Lie group analysis. In problems of MP this feature appeared as 
decisive in creating "RG-algorithm" which has united RG ideology of QFT with a regular 
way of symmetry construction for BVP solutions. Due to this algorithm also there arised 
the concept of "renormgroup symmetry" for solutions of BVP: these symmetries result 
from calculation procedure similar to that used in the modern group analysis. 

At the initial stage [H |2] , application of RG algorithm was mainly limited to problems 
based on DEs though formally this algorithm can be used in any problem, for which it 
is possible to specify a regular way of calculation symmetries for basic equations. Hence, 
transition to such objects, which up to a recent time were not a subject of the group anal- 
ysis, in particular, to the integral and integro-differential equations, essentially expands 
the area of RGS applications. 

In just mentioned cases integral relations form a skeleton of a problem. They, however, 
can appear as some independent objects for application of RGS, constructed for solutions 
of DE. Frequently, of physical interest is not the solution itself in all range of change of 
variables and parameters, but rather its certain integral characteristic, a solution func- 
tional. This characteristic can appear, for example, in result of averaging (integrating) 
over one of independent variables^ or when passing to new integral representation, e.g., 
to Fourier representation. In this case, RG algorithm can be applied not for improving 
of a particular solution with the subsequent calculation of its integral characteristic, but 

1 

^ And representable as the integral J n{x, i)) d cos 19 of the solution of kinetic equation n{x, 19) . 



^see, for example, previous footnote. 
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directly for improvement of the solution functional for an approximate solution. These 
motives gave a stimulus for expanding of the RG algorithm to models with non-local 
(integral) equations. 

The contents of the paper is structured as follows. In section 2 an introductory example 
to RGS algorithm in mathematical physics is presented, illustrated by a solution of a 
simple BVP. In section 3, generalization of RG algorithm, developed in application to 
BVP for DE jH] (and reviewed in j4j and |22j), is formulated for models with non-local 
equations. In section 4, the review of recent results received on the basis of the modified 
RG algorithm is presented, and also efficiency of a method in application to some already 
known solutions is shown. Some generalities, uniting concepts of RG symmetry, the 
functional self-similarity, are considered in Appendix. 



2 The introductory example to the RGS algorithm 

Generally, the RG can be defined as a continuous one-parameter group of specific trans- 
formations of a partial solution (or solution characteristic) of a problem, a solution that 
is fixed by boundary conditions. The RG transformation involves boundary condition 
parameters and corresponds to some change in a way of imposing this condition. 
For illustration, consider transformations T^, 

x' = f{x,a), f{x,ao)=x\ i = l,...,n, (1) 

depending on a real parameter a, where x G M". A set G of these transformations 
form a one-parameter local group if the functions f\x,a) satisfy the composition rule 

f{f{x,a),b) = f{x,(j){a,b)), (j){a,ao) = a , (j){ao,b) = b , (2) 

that can be transformed to the simplest form with 0(a, b) = a + b and with oq = in (^. 

For a given solution of some physical problem renormgroup transformations in the sim- 
plest case are defined as transformations of Q type, i.e. as simultaneous one-parameter 
group transformations Rt of two variables, say x and g, 

x^x' = x/t, g^g' = g{t,g), (3) 

the first being the scaling of a coordinate x (or reference point) and the second - a more 
complicated functional transformation of the solution characteristic g. Hence, the RG 
transformation corresponds to a change in the parametrization for the same solution, 
while the equation Q for the function g has the form 

9{x,g) =g{x/t,g{t,g)) , g{l,g)=g, (4) 

and guarantees the group property Rrt = RrRt fulfillment for transformations Q- These 
are just the RG functional equations and transformations for a massless QFT model with 
one coupling constant !H]- In that case x = Q'^/fJ^'^ is the ratio of a four-momentum Q 
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squared to a "normalization" momentum /i squared and g is a coupling constant, while g 
is the so-called effective coupling. 

Geometrically, transformations (P) mean that any point x G M" is carried by this 
transformations into the point x whose locus is a continuous curve passing through x, 
known as of a path curve of the group G. The group property means that any point of 
a path curve is carried by G into points of the same curve. The locus of the images Ta{x) 
is also termed the G-orbit of the point x. The correspondence between transformations 

and means that for RG transformations a curve in the {x,g} plane that defines 
the solution of a physical problem is the path curve of the renormgroup Rt- In other 
words, the solution of the problem coincides with the i?i-orbit of a boundary manifold - 
the point {x = xo,g = go}. Upon the RG transformations the reference (boundary) point 
{^Oydo} is shifted to some other value {xi,gi}, while the solution remains unaltered, i.e. 
the solution curve is the invariant manifold of the group Rt (like the invariant charge in 
QFT HH). 

Hence, the general problem of searching for the RG transformations may be reformu- 
lated as follows: the solution of the physical problem should coincide with the orbit of 
the renormalization group. 

In mathematical physics a solution of a physical problem usually appears as a solution 
of some EVP. Then the corresponding RG transformations may be obtained from the 
symmetry group related to this EVP, provided the boundary condition is also involved in 
group transformations. The key point here is the fact that the corresponding symmetry 
group are calculated using the regular algorithms of modern group analysis, provided the 
basic mathematical model is formulated in terms of differential (or integro-differential) 
equations. 

Let this model be given by a system of k—th order differential equations, identified 
with its frame, 

F^(x,M,M(i), . . . = , a = l,...,s. (5) 

In the paper we use the terminology of differential algebra and notations for variables 
accepted in the group analysis |24j: 



X = {x'}, U = {u"}, = {<}, M(2) = . . . , (6) 

where a = 1, . . . , m; i, j, . . . = 1, . . . ,n. Variables x and u are referred to as independent 
variables and differential variables, respectively, having consecutive derivatives U(2), 
. . . etc. Differential variables are related by a system of equations 

< = A («") , «° = «) = A («") , . . . (7) 
via the operator of the total differentiation 

Locally analytical function of variables (jH)), for example, the function F{x, u, ^(i), . . . , U(k)) 
with the highest order derivative k refers to as differential function of the k-th order, and 
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a set of all such functions with any values of k forms the space of differential functions 
^[XjM]. Any function F G gives rise to a differential manifold [F], determined by 

an infinite system of equations 

[F]: F = 0, DiF = , DiDjF = , . . . . (9) 

The manifold [F] is called the frame of the k-th order partial differential equation 

According with the definition a system of s-th order differential equations is said to 
be invariant under a group G if the frame of the system is an invariant manifold for the 
extension of the group G to the s-th order derivatives |25| p. 209]. When utilizing an 
infinitesimal group generator 

X = eU.«+r7-5„., Cv'^eA, (11) 

with coordinates ^* , 77", which are functions of group variables {x^,u°'}, this definition 
leads to the invariance criterion in the following form 

=0, (r = l,...,s, (12) 

where X(^k) denotes X extended to all derivatives, involved in F^r and the symbol \ [f] means 
evaluated on the frame Q. Solving a system of linear homogeneous partial differential 
equations (known as the determining equations) for coordinates gives a set of in- 

finitesimal operators (jllll (or group generators) which correspond to the admitted vector 
field of the symmetry group G and form a Lie algebra L. 
Let the Lie group G with the generator 

X = edt + Cd, + vdy, (13) 

be defined for the system of the first order partial differential equations 

Vt = F{t,x,y,y^) . (14) 

The typical BVP for (jl4|) is the Cauchy problem with the boundary manifold defined by 

t = 0, y = ^{x). (15) 

The solution of this Cauchy problem is the G-invariant solution, iff for any generator (fT^ 
the function ip obeys the Equation ^26,, §29]: 

77(0, X, ^) - r(0, X, V^)9,V^ - e*(0, X, V^)F(0, X, V^, d^^) = . (16) 

The solution of the Cauchy problem (fT^ - (fT3j) coincides with orbit of the group G and the 
boundary manifold is partially invariant manifold of the group G with the defect 6 = 1. 
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This example gives an instructive idea of constructing generators of renormgroup sym- 
metries. The milestones here are: a) considering the BVP in the extended space of group 
variables that involve parameters of boundary conditions in group transformations, b) 
calculating the admitted group using the infinitesimal approach, c) checking the invari- 
ance condition akin to ()16p with the goal to find the symmetry group with the orbit which 
coincides with the BVP solution and d) utilizing the RG symmetry to find the improved 
(renormalized) solution of the BVP. 

The full algorithm was described in details in our previous publications [HI Ej and will 
be also touched upon in the next section while here we will only give a general grasp at 
the problem using a trivial example, the BVP for the Hopf equation 

v, + vv^ = 0, v{0,x) = eU{x), (17) 

where U is the invertible function of x. Introducing u = ev we insert the boundary 
amplitude directly in the input equation 

Uz + euux = , u{0,x) = U{x). (18) 

For small values of ez <^ 1, i.e. near the boundary, z — > 0, or for small amplitude at 
the boundary, e — 0, a perturbation theory solution to (jTH)) has the form of a truncated 
power series in ez, 

u = U- ezUUx + O {(ezf) . (19) 

It is obvious that this solution is invalid for large distances from the boundary, when 
ezUx — 1. The renormgroup symmetry gives the root to improving the perturbation 
theory result and restore the correct structure of the BVP solution. 

With the goal to obtain this symmetry we extend the list of variables, involved in 
group transformations, adding the parameter e to the list of independent variables. Then 
we calculate the admitted symmetry group G, with the generator 

X = Cdz + Cd. + Cde + vdu, (20) 

using the classical Lie calculational algorithm (see, e.g. [22]) that employs the infinitesimal 
criterion ()12p. The solution of determining equations gives coordinates of the generator 

(EOD, 

^' = ij\ r = ewT/;^ + ^ x(V^3 + ^^) , r = e^\ r] = utp^ (21) 

where tp^ , i = 2,3, 4, are arbitrary functions of e, u, x — euz and ip'^ is an arbitrary function 
of all group variables. These formulas define an infinite-dimensional Lie algebra with four 
generators^ 

Xi = ^P' (a, + eudx) , X2 = iJ^dx , = V'' (xd, + udu) ,X^ = ^^ (e9, + x9,) . (22) 

•^In case when the amphtude e is not involved in transformations we have only three generators (see, 
e.g. El p.222]). 
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Suppose we are given a particular solution of the BVP (fTHj) . u — W{z,x,e) = 0, which 
defines an invariant manifold of the group (j2(Jj) . (j2H) . The corresponding invariance con- 
dition, evaluated on the frame (fTH|) . looks similar to (fT2|) . 

{W - xW^)^^ - W^^"^ - {eW, + xW^)^^ = . (23) 

This equation is valid for all z, hence it remains valid for z = 0, when W is replaced by 
U{x). In this limit, z ^ 0, the condition ()23p gives a relationship between ip\ i = 2,3, 4, 
(no restrictions are imposed on ip^), that can be easily prolonged on z ^ 0, 

i^' = -X{^' + + (w/f/x)V'' , X = x-euz, (24) 

where the derivative should be expressed either in terms of x or ti in account of the 
boundary condition. Inserting ()24|) in !^21\i we get the group of a smaller dimension with 
the generators 

Ri = (5. + eud^) , R2 = ui;^ [{ez + d, + 9J , = ez^^ {zud, + d,) . (25) 

The above procedure, that transforms to we refer to as the restriction of the 
group (j2(jp on a particular solution. 

The solution of the BVP defines a manifold, which appears to be invariant for any 
generator Ri just by a method of construction, hence ()25|) defines the desired RG symme- 
tries. This means that the solution of the BVP can be constructed using any generator of 
the RG algebra (jSHI), say, the generator R3. Without loss of generality we choose eip'^ = 1 
and obtain the following finite RG transformations (a is a group parameter) 

x' = X + azu , e' = e + a , z' = z , u' = u . (26) 

Here z and u are invariants of RG transformations, while transformations of e and x are 
translations, which in case of x also depend upon z and u. For e = we have, in view of 
(fT^ . X = H{u), where H{u) is defined as the function inverse to U{x). Then eliminating 
a, z and u from ()26p and omitting primes over the variables, we get the desired solution 
of the BVP (|18|) in the non-explicit form 

x-ezu = H{u) . (27) 

This in fact is the improved PT solution (jl9p . which is valid not only for small e <^ 1, 
provided the dependence (j77j) can be resolved in a unique way. 

The peculiarity of the procedure of constructing RG symmetries is the multi-choice 
first step that depends on the way in which boundary conditions are formulated and 
the form in which the admitted symmetry group is calculated. For example, instead of 
calculating the Lie point symmetry group we may consider the Lie-Backlund symmetries 
[77] with the canonical generator R = kO^, where k depends not only on z, x, e, u, but on 
higher order derivatives of u as well. We may look for k in the form of a power series in 
e, and the invariance condition (j^^ is formulated as vanishing of k at z — 0. Depending 
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on the choice of the zero-order term representation we get either infinite or a truncated 
power series for hnear in e form, 



R = ndu, K = l - f , - ezUx . (28) 

This RG generator ()28j) is equivalent to the Lie point generator R2 from ()25j) . thus giving 
the same resuh. 

Another possibihty of calculating RG symmetries for BVP (fTH|) is offered by taking 
into account some additional differential constraints which are consistent with the bound- 
ary conditions and input equations. For example, if the boundary condition in (fTHj) is 
linear in its argument, U{x) = —x, the differential constraint may be taken as Uxx = 0; 
this equality reflects the invariance of the original equation with respect to the second- 
order Lie-Backlund symmetry group. Calculating the Lie point symmetry group for the 
joint system of this constraint and the Hopf equation gives another way to flnding RG 
symmetries for the BVP fjl8|) . 

The above example demonstrates the milestones of the RGS method in mathematical 
physics. The details of the general approach will be discussed in the next section. Here 
we just point to the fact that in order to construct RGS we employ the symmetry group 
calculated in a regular way using the modern group analysis technique. This group is 
considered in the extended space of variables that includes the parameters and boundary 
data involved in group transformations. The invariance of the perturbative solution is used 
to flnd the particular RGS generator that leaves the particular BVP solution unaltered. 
Then the utilization of the flnite group transformations restore the desired structure of 
the BVP solution. 



3 The scheme of RG algorithm for non-local prob- 
lems 

A procedure of construction and use of RGS with reference to BVP for DE is described in 
details in our previous works, for example, in reviews jSJ and illustrated in the previous 
section. Nevertheless, for connectivity of a statement, the basic stages of the scheme are 
briefly depicted below with accent on those changes, which it is necessary to introduce in 
RG algorithm in order to make it applicable to non-local problems. 

The starting point in these problems is the system of z/ > 1 integro-differential (includ- 
ing differential and integral) equations for functions u = a = 1, . . . , m of variables 
X = {x*}, i = 1, . . . ,n, 

[E]: E,iuix)) = 0, (29) 

with the non-local terms depending upon integrals of these functions, and supplemented by 
appropriate boundary or initial conditions. We assume, that we know some approximate 
solution, = f/", for example, represented by a truncated power series in powers of a 
small parameter or in powers of a small deviation from a boundary of an area, where the 
solution is known. 



8 



Perturbative 
BVP solution 



Basic 
model 



Basic manifold HAd 



II 



Generators of group 



III 



RG invariant solution 







RGS generators 










IV 





Figure 1: General scheme for realization of RG algorithm. 



Then, conditionally, the scheme of realization of RG algorithm can be expressed as a 
sequence of four steps, submitted by the scheme on the Figure HJ 

(I) construction of basic manifold, 

(II) calculation of the admitted symmetry group and 

(III) its restriction on a particular BVP solution leading to revealing of RGS, and also 

(IV) searching for an analytical solution that is adequate to this RGS. 



3.1 Construction of RG manifold 

The initial problem is the construction of RGS and appropriate transformations that also 
touch upon parameters of particular solutions. Therefore the purpose of a first step (I) 
consists of involving in group transformations in this or that way the parameters entering 
both in the equations of a problem, and in boundary conditions on which this partial 
solution depends. This purpose is achieved by construction of a special manifold TZAi 
which we believe is given as a system of s DE of the fc-th order and q non-local relations, 

F„{z,u,U(i), . . . ,U(k)) = , a=l,...,s, (30) 

F^{z,u,ui^i),...,u^r),J{u))=0, a = l + s,...,q + s. (31) 
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Non-local variables J{u) included in (jHT|) are introduced by integrated operations 



A presence of relations (j3ip characterizes the basic difference of 71A4 for non-local prob- 
lems from the case of BVP for DE, for which the manifold TZAi is differential. 

Let's note that generally TZA4 does not coincide with a system of initial equations. 
Only on occasion, which it will be told about below, and at the additional clauses expand- 
ing and specifying the list of variables and the parameters entering in RG transformations, 
it is possible to establish conformity between and (IHU jl - lj^ . However it is not pos- 
sible to execute a first step of the algorithm for any boundary problems simply treating 
[E] as TZAi, and the concrete form of a realization of the first step depends both on a 
form of the initial equations ()29p . and on the form in which boundary conditions are 
presented. Formulated earlier for models with DE the general approach to construction 
of TZAi remains also valid for non-local problems. In this sense several specified earlier 
(see, for example, 0) and illustrated in section I ways of constructing this manifold are 
possible. These ways do not exhaust all opportunities to construction of TZAi, rather 
they emphasize a variety of approaches to realization of a first step depending on a char- 
acter of a problem under consideration. A choice of a concrete realization more often is 
dictated both by a form of the Eq. and boundary conditions to them, and a type of 
the approximate solution. Such multi-variant situation is inherent only in the first step 
of the algorithm and is aimed on covering the widest variety of the problems, analyzed 
by this method. Already the following step of the scheme is carried out in frameworks of 
the well-developed group-theoretical methods. 

3.2 Calculation of the transformation group 

The next step (II) consists in calculation the most wide admitted symmetry group Q 
for Eqs. ()30 |) -()31 |) . Here the essential change of RG algorithm is required in comparison 
with its realization for differential manifold TZAi. Really, in application to TZAi, defined 
only by a system of DE (j3U|) . the question was about a local group of transformations in 
space of differential functions A, at which system (jHU)) remains unchanged. 

At transition to manifold TZAi which is set by the system of Eqs. fj30p -()3H ) . classical 
Lie algorithm, using the infinitesimal approach, appears inapplicable. The basic obstacle 
here is that TZAi in this case is not determined locally in space of differential functions, 
therefore the basic advantage of Lie computational algorithm, namely, representation of 
DetEq as the over-determined system of equations is not realized here. Also in frameworks 
of the classical group analysis the procedure of prolongation of the group operator of point 
transformations on non-local variables is not defined. Probable ways to overcome these 
complexities while performing the second step of RG algorithm are specified below. 

At modification of the RG algorithm we lean on the direct method of calculation 
of symmetries which was advanced in I2H|^I2H] and used for finding symmetries for the 




(32) 
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Boltzmann kinetic equation, the equations of motion of viscous-elastic media and Vlasov- 
Maxwell equations in the kinetic theory of plasma. This method is based on a generaliza- 
tion of a group of symmetry, the so-called Lie-Backlund symmetry group (terms "higher" 
or "generalized" symmetry are also in use), defined by the generator of the form (jlip 
prolonged on all higher-order derivatives, 

X = Cd,. + tdu^ + Qd^ + c^,9„. + . . . , 

c = A(x")+e<,., c^., = A.A.(x")+ew^,^,,, x° = r/"-e<, 

with coordinates ^'^{\z,u\), rj°'{\z,u\)^ Ci*([-2? w]) . . ., being differential functions from the 
space A. A set of all Lie-Backlund operators forms an infinite-dimensional Lie algebra Lg, 
and any operator of a form X^, = C,^Di is the Lie-Backlund operator for any differential 
function ^''{[z,u]); the set of operators X^ forms an ideal in Lg. This property allows 
to introduce a notion of equivalence of two Lie-Backlund operators Xi , X2 G if 
Xi — X2 G (written as Xi ~ X2). In particular, any Lie-Backlund operator X G is 
equivalent to the operator (jHHI) with = 0, 

X ~ r = X - f A = x°5n- , = r]" - Cut ■ (34) 



The operator Y is known as the canonical representation of X, and in notation we 
imply the prolongation of action of the operator on all higher-order derivatives according 
to formulas (j33j) . It is essential, that in the group of infinitesimal transformations Q with 
operator and the parameter a only dependent variables are changed. 



m'" = + ax" + O(a') , z" = z' , (35) 

while independent variables z^ remain unchanged. This property has allowed to formulate 
the concept of symmetry groups of IDE of form ()3H) as a local group of transformations Q 
with operator ()34j). at which the form of the function Fcr remains unchanged for any value 
of a group parameter a. Differentiation of the appropriate invariance condition, which 
has been written down for the function F^- dependent on the transformed dependent 
variable u'", with respect to the group parameter a and transition to a limit a ^ 
gives determining equations (DetEq). In difference from a case of basic DE, these DetEq 
generally also are non-local. 

With the help of the canonical operator Y, the invariance criterion for Eq. (j31|) with 
respect to the admitted group can be written down in an infinitesimal form 



YF. 



0, a = l + s,...,q + s, where Y = [ dz >c(z) ^ . (36) 

J ou {z) 



Meaning generalization of action of the canonical group operator not only on differential 
functions, but also on functional, here in definition of Y variational differentiation [21] 
is used. For integral functional (jH^ a derivative of 6J/6u (z) with respect to a function 
u is defined via the principal (linear) part of an increment of a functional as a limit (if it 
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exists) (see PT|1 



5J\v\ , J \u + 5Uf\ — J \u\ 
5u [z) e^o dr (r) 



Here an infinitesimal variation 5u^{z) > is a continuously differentiable function given 
on fixed interval A = [ri, T2\ which differ from zero only in e- vicinity of a point z, and 
the norm ^ at e ^0. 

It is checked up by direct calculation that the action of the operator Y on any dif- 
ferential function and its derivatives, for example u^Uz-, ■ ■ ■ gives a usual result: Yu = 
X, Yuz = Dz{>c) etc. Hence, if = is a usual DE, formulas result to local DetEq, 
and in a case when = has a form of a system of IDE, formulas can be considered 
as non-local DetEq, dependent both upon local and non-local variables. 

Treating local and non-local variables in DetEq as independent variables allows to 
divide these equations into local and non-local. A procedure of a solution of local DetEq 
is carried out in a standard way, using Lie algorithm based on splitting of a system 
of overdetermined equations with respect to local variables and their derivatives. In 
result expressions for coordinates of the group operator are found, determining so-called 
group of intermediate symmetries ^29j which are used further at the analysis of non- 
local DetEq. Procedure of the solution of non-local DetEq is carried out similarly, by 
substituting coordinates of group operator of intermediate symmetry found in non-local 
DetEq and splitting them with the help of variational differentiation. Hence, construction 
of symmetries for non-local equations also becomes an algorithmic procedure, which can 
be presented as a sequence of the following operations: 

a) definition of a set of local group variables, 

b) construction of DetEq on a basis of infinitesimal invariance criterion, 
which uses a generalization of the definition of the canonical operator, 

c) separation of DetEq into local and non-local, 

d) solution of local DetEq with use of the standard Lie algorithm, 

e) solution non-local DetEq with the help of the operation of the variational differen- 
tiation. 

These operations are generalization of the second step of algorithm on a case, when TZM. 
is the integral or the integro-differential manifold. 

In summary we describe an operation of prolongation of a symmetry group on non-local 
variables. To execute standard (in the classical group analysis) operation of prolongation 
of the operator of Lie point group on the non-local variable defined, for example, by 
integral relationship (jH^ . first write down this operator in the canonical form, F, and 
then formally prolong it on non-local variable J 

Y + x^dj ^ xdu + >t^dj . (37) 
^Definition of a variational derivative for arbitrary functional see in [221 
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The integral relation between x and is obtained by applying operator ()37|) to Eq. (jH^ 
which was introduced as a definition of the variable J. Substituting explicit expressions 
for the coordinate h of the operator Y and calculating the resulting integrals we obtain 
the required coordinate x'^ of the prolonged operator, 

xiz) dz= f ^'^/''/f^^ x(z) dz dz' = I Mz) dz . (38) 
J (^u{z) J 

Here for brevity only one argument of a generator's coordinate is specified, namely the 
argument upon which the integration is fulfilled. 

The actions described in this paragraph resulting in operators of the admitted group 
in non-canonical (jHHj) or in canonical ()34|) representation, make essence of the second step 
of RG algorithm. 




3.3 Restriction a group on a solution 

The group found on the second step, Q, which is determined by operators ()33|) and 
is generally wider, than the renormalization group of interest. The last is related to the 
concrete particular solution of a boundary problem, hence with the goal to get RGS it 
is necessary to make the third step (III), consisting in restriction of the group Q on 
a manifold, set by this particular solution. From the mathematical point of view, this 
procedure consists in checking vanishing conditions for a linear combination of coordinates 
x" of the canonical operator equivalent to ()33|1 . on some particular solution U°'{z) of a 
boundary problem 



i j ) 



= . (39) 



The form of the condition set by relation (jH^ . is common for any solution of the EVP, 
but a way of realization of the restriction procedure of a group in each separate case 
is different. Usually a particular approximate solution for a concrete EVP is used. On 
the general scheme it is specified as the dotted arrow connecting "initial object" — the 
approximate solution for particular EVP - with that object which arises in a result of the 
third step. 

At calculation of combination (jH^ on a concrete particular solution U°'{z) it is trans- 
formed from a system of DE for group invariants to algebraic relationships. We shall 
note two consequences of the specified actions. First, the procedure of restriction results 
in a set of relations between various and thus "links" coordinates of various group 
operators Xj, admitted by TZJ^ finn|) - (jHT|) . Second, it eliminates (in part or completely) 
an arbitrariness which can arise in values of coordinates 77" in a case of the infinite 
group Q. 

As a rule, the procedure of restriction of the group Q reduces its dimension. At that 
the general element (j33p of this new group TZQ after performance of this procedure is 
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represented by a linear combination of new generators Ri with coordinates 77", 



X ^ R = J2 B^Rj ' = ijdx^ + Vjdw 



(40) 



and arbitrary constants BK 

The set of operators Rj, each containing the required solution of a problem in the 
invariant manifold, defines a group of transformations, which by analogy with RG for 
models with DE we also refer to as renormgroup. 

The statement of the third step of the algorithm given in this section finishes the 
description of a procedure of construction of RGS. In the following section it is shown, 
how RGS are used for achievement of an ultimate goal of the RG algorithm, namely 
improvement of an approximate solution. 

3.4 Construction of RG-invariant solution 

Three steps described above completely define the regular algorithm of construction of 
RGS, but one more, the fourth (IV), final step is necessary. This step consists in use of 
RGS operators for finding analytical expressions for new, improved (in comparison with 
initial) BVP solutions. 

From the mathematical point of view realization of this step consists in use of renorm- 
group invariance conditions^ which are set by a joint system of the equations (jHUjl - pH) and 
vanishing conditions for a linear combination of coordinates of the canonical operator 
equivalent to (pUj). 



Necessity of use TZAi while constructing BVP solution is marked on Figure by the 
dotted arrow connecting these two objects. 

It is clear, that the form of ()4H) is close to ()39p. However, contrary to a previous 
step, differential variables uf in ()41|) are not replaced with the approximate expressions 
for BVP solutions U{z) but are treated as usual dependent variables. 

In the most widespread case, when the renormgroup appears as a one parametric Lie 
point group, the invariance conditions yield first order partial differential equations. Solu- 
tions of the connected characteristic equations give group invariants (similar to invariant 
charges in QFT) through which the required solution of BVP is expressed. Generally for 
arbitrary RGS the RG invariance conditions, written down for any BVP, do not represent 
characteristic equations for the Lie point group operator. They can have more complex 
form, for example, being represented as a combination of partial DE and ordinary DE of 
a high order. However the common approach to construction of BVP solution as the RG 
invariant solution remains in force. 

■''Here we should point on existing the analogy between this concept and concept of functional self- 
similarity (see also discussion in the Appendix inj . 




(41) 
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The description of the fourth step finishes the description of the regular algorithm of 
construction of RGS for models with IDE. It should be noted that two last steps are in 
the same root as for models with DE. The following section contains a number of concrete 
examples showing the ability of RGS algorithm. 

4 Construction of RGS for integral models 

4.1 The example with functionals of the solution of Hopf equa- 



Let's proceed with a simple illustrative example, that we have discussed earlier at the 
beginning of the paper, i.e. an initial problem for Hopf equation (jl7|) . We have shown 
that the solution of the EVP (fTTj) can be constructed using any generator of the KG 
algebra (j^ . Let we are interested not in the solution in all space, but only in some 
solution characteristic at a specific point, say, a value of its first derivative at a point 
X = 0, which formally can be introduced by a linear functional of u, 



The dependence of upon z can be easily restored using the prolongation of the linear 
combination of RG generators ()25|1 on the solution functional (02)) • We use again the 
last generator from the list ()25p in its simplest form with eip'^ = 1. Then we write down 
this generator in the canonical form and calculate its prolongation using formulas ()37p 
and (|nHj). Restricting the RG operator obtained after such prolongation to the space of 
group variables {2;,e, we get the RG generator for the solution functional (ji^ . For 
concreteness we choose U = —x and get the following generator 



The initial condition for the m° at 2; = is known, vP^{z = 0) = —1, hence the use 
of the invariant of this generator, J*^ = ez — = 1 restores the desired dependence 
M° = —1/ (1 — ez), which is valid from the point 2; = up to the singularity point 
Zsing = l/c- Wc uotc that this result is obtained without construction of BVP solution, 
using only the appropriate RGS. On the first sight the considered methodical example 
and the construction carried out looks clumsily enough and it is much easier to proceed 
from the trivial solution ()27|) . But in more complex situations the explicit form of the 
solution is frequently unknown, whereas it is possible to construct RGS. In two subsequent 
sections of the paper, 14.21 and 14.31 we give examples that serve as an illustration of this 
statement. 



tion 




(42) 



(43) 
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4.2 Nonlinear optics: development of a singularity on the 
laser beam axis 

We proceed our illustrations with analyzing the BVP for a system of two first order partial 
DE, which include Hopf equation as a particular case and are widely used in gas dynamics, 
optics and plasma physics: 

Vz + vvx-alx = 0, h + vh + Ivx = vlv/x . (44) 

For concreteness below the terminology of nonlinear optics is used where (j44j) are known as 
nonlinear geometrical optics equations (see, for example, the review [23) and are utilized 
to describe the evolution of a laser beam in a nonlinear medium. In this case I{z,x) 
stands for the beam intensity and f (z, x) is the derivative of a beam eikonal with respect 
to transverse coordinate, a is a factor of nonlinear refraction, z and x are coordinates in 
a direction of a beam propagation and in a perpendicular direction, respectively; u = 1 
and = for cylindrical and for plane case, respectively. 

The nonlinear medium occupies a half-space z > and boundary conditions for equa- 
tions are set at z = 0, 

v{0,x) = I{0,x)=I{x). (45) 

The specific choice of the zero value of an eikonal derivative on the boundary corresponds 
to a collimated beam with the distribution of the beam intensity upon a transverse coor- 
dinate X, characterized by the function T{x) with the maximal value equal to unity. 

Various analytical approaches (see the papers f3^, TT]) do not provide us with a 
universal method for finding solutions to BVP (j44|l - (j45|l . suitable for any geometry of a 
problem and any boundary data, while the use of RG algorithm offers new opportunities 
(see, for example, [3H1 II])- In particular, it allows to prolong the solution, known in 
a small vicinity of the boundary z ~ of a nonlinear medium, up to a vicinity of a 
solution singularity, which occurrence represents the most attracting physical effect. The 
singularity is formed on a beam axis hence the behaviour of I{x) and f (x) at a; — *• is of 
particular interest. 

The surprising thing is that this behaviour can be understood without knowledge of 
the complete solution via the application of RG algorithm to two functionals of the BVP 
solution namely, to the intensity of a laser beam I'^{z) = l{z,0) and to the 

second derivative of the eikonal W^{z) = Vx{z,0), calculated on an axis of a beam and 
related to this solution by formal relationships 

/»W^/d. w»w^/d..(.K(...). (46) 

Boundary conditions for these functionals with the account of (|45p are given as 

/°(0) = 1, W\0) = 0. (47) 

In spite of the fact that these conditions do not contain data on the dependence of a beam 
intensity upon the coordinate x, such information is included in the RGS operator which 
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explicit form is defined by the profile of the beam intensity I{x) ai z = 0. We present two 
examples, corresponding to cylindrical and plane ("slit") laser beams with various X{x). 

For a cylindrical beam {u = 1) with the parabolic intensity distribution I{x) = 1 — x^, 
the RGS generator has the form 

j^par = (1 _ 2az'^) d, - 2azxd^ -2a{x- vz) + Aalzdi . (48) 

To determine the dependence of /° and on the coordinate z we prolong ()48|) on non- 
local variables (solution functionals) /° and that gives the following generator in the 
reduced space of variables {z, J°, 

i?4 = (1 - 2az^) d, + Aal^zdjo - 2a{l - 2zW°)dwo . (49) 

The use of two invariants of generator Ji = {l — 2az'^)I^ and J2 = W^{l — 2az'^)+2az 
with evident equalities Ji = 1 and J2 = 0, which follow at the account of boundary 
conditions (|T7j). immediately gives expressions 

' W' = --^. (50) 



1 - 2az^ ' 1 - 2az'^ 

These formulas describe spatial dependence of variables I^{z) and W^{z), starting from a 
boundary of a nonlinear medium z = up to the point Zsing = l/V2a, where the solution 
singularity occur, i.e. where the beam intensity and the eikonal derivative turns to infinity; 
beyond this point there is an area of rays intersection, where equations can not be 
applied. Expressions (j5(J|) also follow from formulas, obtained earlier [321 without use of 
RG algorithm. However, the RGS algorithm here presents an elegant way of obtaining 
these formulas without calculating the complete solution to EVP. 

Curves of typical dependencies of variables I'^{z) and W^{z) upon the dimensionless 
coordinate z/zsmg at a = 0.1 are given on Figure |21 The change of the parameter a does 
not change a type of the curve for the intensity /°, whilst values of on the right panel 
vary proportionally to y/a. Block curves correspond to formulas ()5()|1 . i.e. to parabolic 
distribution of intensity of a cylindrical beam at the medium boundary; dotted curves 
refer to plane geometry of the beam, considered below. 

The procedure of prolongation of the operator on non-local variables uses a canonical 
form of RG generators (Lie-Backlund operators) and is suitable also in that case when 
this generator is given by a higher-order Lie-Backlund symmetry. Such case is realized 
for a plane laser beam with "soliton" profile of the intensity distribution at the boundary. 
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0.2 0.6 1 ^sing 0.2 0.6 1 ^sing 



Figure 2: Dependencies of the intensity of a laser beam (at the left) and the second 
derivative of its eikonal (on the right) on the beam axis x = at various distance zj z^mg 
from the boundary of a nonlinear medium 2r = 0, plotted with the use of formulas for the 
cylindrical (block curves) and plane (dotted curves) geometry. 



T(x) = cosh ^(x), when the appropriate RGS generator has rather cumbersome form 

~ ' ' 2 



R 



sol 



{Ivl + allY 



^ [Ivl - "4^) (t^^ + 4a(l - /)) + ^avll^v^ v^^ 



+ {2av [all - Ivl) + c^^xh (t;' + 4a(l -/)))(/,,- ^ 



v{l - tv^) - atl^ 



+ 



- [Ivl - all) [v"" + 4a(l - J)) + 4avlvj, ) ( 4. - ^ 



1 



[2v [all - Ivl) + vJx [v^ + 4«(1 - /))] Iv^, + —{Ivl + [4aZ 



+ {hv - 21v:,f] ] - /(2 - tv^) + tv^Y^ . 



(51) 

Prolongation of (jSlj) on non-local variables (j46|) gives the more simple operator in space 
of functionals 



i?5 = K^'^djo + x^°dwo 



(^4-5/°-z/° + 2(/°-l)^j9,o 



JO jO 



JO 



- 2(/° - 1) 



2(4^ 



\2 



(J0)2 (J0)2 (J0)3 



(52) 



While obtaining this formula we used the relation between derivatives of functions / and 
Vx with respect to spatial variables on the beam axis (at x = 0), which follows from the 
initial equations, 

rO\ 3 jo jo 



Vxxx^^y 0) 



aP 



P 



+ 10 - 



P 



0\2 



P 1 

^x(2,0) = -^, lxxiz,0) = - 



JO 



It 



2{ ^] - ^ 



P 



(53) 
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The beam intensity and its second eikonal derivative on an axis are defined from RG 
invariance condition (j41|) . wliicli is equivalent to vanisliing of coordinates and 
of the operator (j52j) . This condition gives two ODE of the second and the third order 
respectively. Solving at first ODE of the second order with the initial conditions (jTTj) and 
the additional condition on the first derivative, (/°/a/J° — 1)1^^0 = 2^/^ which follows 
from the relation (jSHj) at 2 = 0, we obtain in the implicit form the law of variation of 
and (compare with formulas (|5U|) for a parabolic beam), 



These formulas are valid from the boundary of a nonlinear medium 2; = up to a point 
where the solution singularity occurs. The coordinate of the solution singularity is found 
in view of the fact that the derivative of turns to infinity at this point, that gives 



z 



smg 



1/21/0, and the value of the intensity in this point is equal to two. The solution 
of the remaining ODE of the third order gives the same result. Dependencies of and 
upon the dimensionless coordinate z/zsmg at a = 0.1 are plotted on Figures El by 
dotted curves. Without prolongation of RGS on non-local variables the result ()54|1 with 
the use of RG algorithm was obtained earlier in |3] , though in a more complicated way. 

Summarizing this paragraph, we note, that universality of a procedures of prolonga- 
tion of RG generators presented either as point group operators, or Lie-Backlund group 
operators has allowed to describe from uniform positions an occurrence of a singularity 
of the BVP solution to using for this purpose the reduced description in terms 

of solution functionals. 



4.3 RGS for solution functionals of plasma kinetic equations 

In the paragraph we continue with demonstrations of potentialities of RGS algorithm 
for nonlocal models. In contrast to the previous paragraphs, where we described the 
prolongation of the RGS generators on solution functionals for differential models, here we 
consider the case when integral relations form a skeleton of a problem. A vivid example 
is the model that is used in the plasma kinetic theory and define the evolution of a 
collisionless inhomogeneous plasma. 

The macroscopic state of plasma particles is governed by distribution functions /" 
(for each species of plasma particles with mass and charge Cq,), that dependent on 
time t, a coordinate x of a particle, and its velocity v (for simplicity we consider the one- 
dimensional plane geometry). Evolution of distribution functions is described by Vlasov 
kinetic equations [10] . 

/r + vf: + {eJm^)E{t, x)f: = (55) 
supplemented by Poisson and Maxwell equations for the electric field E, 

E,-AttJ2 [ dve^r = , Et + AnY, I dwe^r = • (56) 
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The joint system of equations and the first equation in (jSUj) is often referred to as 
Vlasov-Poisson (VP) equations. We are interested in a solution to the Cauchy problem 
to equations (j55|) with the initial conditions that correspond to the electron and ion 
distribution functions specified at t = 

r\,^, = f^ix,v). (57) 

The admitted symmetry group for ()55|) . calculated by the method prescribed in 

section \'A.2\ is given by time and space translations, Galilean boosts and two generators 
of dilations 30^. VP equations seem to be the simplest one-dimensional mathematical 
model, which is commonly used to describe the evolution of inhomogeneous plasma, e.g., 
the expansion of a plasma slab. Even so analytical methods fail to create the spatially 
symmetric solution of ()55 |) -()56 |) for the distribution functions with initial zero mean ve- 
locity. Thus, with the goal to find physically reasonable solution we are forced to simplify 
the basic system of VP equations. 

One possible way to simplify the system (j^UI) is to study dynamics of plasma 

expansion in quasi- neutral approximation |1TJ E], suitable for a description of plasma 
fiows with characteristic scale of density variation large as compared with Debye length 
for plasma particles. It means that one can neglect the field terms in Poisson and Maxwell 
equations (fKHjl and consider the total charge and current densities equal to zero. Hence, 
particle distribution functions /"(t,x, f) for the electrons (a = e) and ions (a = 1, 2...) 
obey kinetic equations ()55|) and are assumed to satisfy the non-local quasi-neutrality 
conditions, 

/"dt;5^e„r = 0, /"dt;t;5^e„r = 0, (58) 

a a 

while the electric field E is expressed in terms of moments of distribution functions: 



(59) 



Analytical study of such yet simplified model represents the essential difficulties, but due 
to application of RG algorithm it is possible not only to construct solution at various 
initial particle distribution functions 42 but also to find particles density and energy 
spectra without calculations of distribution functions for particles in an explicit form. 

To construct RGS we consider a set of local and non-local equations as TZAi, 
in which the electric field E{t, x) appears as an unknown function of the coordinate x and 
time t. The Lie group of point transformations admitted by this manifold is calculated 
in a way similar to that used in section 13.21 Here, besides time and space translations, 
the Gahlean boosts and three operators of dilations, there arises a new projective group 
generator jlHj. Precisely this generator enables to construct a class of exact solutions to 
the initial problem that are of interest, as a linear combination of the generator of time 
translations and the projective group generator leaves the approximate PT solution of the 
initial value problem /" = JqIx, v) + 0{t) invariant at t — > 0, i.e. it is the RGS operator 

i?6 = (1 + nH^)dt + nHxd,^ + n^x - vt)d^ . (60) 
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The generator (jUUI) is the only which selects the spatially symmetric initial distribution 
functions with zero mean velocity. The value fl can be treated as the ratio of the ion 
acoustic velocity to the gradient length Lq. 

Group invariants of the RG operator ()60p are particle distribution functions /" and 
combinations J3 = x/yl+TPt^ and J4 = f ^ + Q'^{x — vty. Hence, EVP solutions, i.e. 
distribution functions at any time t 7^ 0, are expressed with the help of these invariants 
in terms of initial distributions (j57j) . 

r = /o"a^"^) , /^"^ = + —$0(^3) . (61) 

Here the dependence of $0 upon the invariant J3 = x/vT+TFt^ is defined by quasi- 
neutral conditions ()58|1 . and the electric field E = —^x is found with the help of a 
potential 

<!>{t,x) = <l>o{J3){l + nH')-' . (62) 

A variety of examples, illustrating these formulas for the plasma slab consisting of different 
groups of hot and cold electrons and ions of various species, is found in • 

Distribution functions ()6H) give exhaustive information on the kinetics of plasma bunch 
expansion. However, for practical applications rough integral characteristics, such as 
partial ion density, n'^{t, x), {q = 1, 2, . . .) and ion energy spectra, dNg/de, might be more 
useful, 

71'^= [ dvnt,x,v), ^ = J- [ dx{nt,x,v) + nt,x,-v)) . (63) 

J de ruqv J 

— oo — oo 

In view of the complex dependence upon the invariant /'-"^ it is not always possible 
to carry out direct integration of a distribution function over velocity in the analytical 
form, therefore here the procedure of prolongation of the operator on solution functionals 
described in section ITT] comes to the aid. The density n'^{t,x) is a linear functional of f^, 
hence we prolong the generator Rq on the solution functional (jUHj) to get the following 
RG generator in the space of variables {t,x,n'^}, 

i?7 = (1 + nh^)dt + nhxdx - nHn'^dn, . (64) 



Two invariants of this generator, namely J3 and J4 = n'^y/l+TWt^ are related for arbitrary 
t 7^ via their initial values: J3|t=o = ^'y '^4|f=o ~ -^gi^')- Therefore, we get formulas 
that characterize spatial-temporal distribution of the density of ions of a given species in 
terms of the initial density distribution 



n 



1 



, Nn (x) , K= , X = , ■ (65) 
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The general form of dNg/de is rather comphcated but its asymptotic behavior at fit — »• oo 
is described by the same function Afg, 



dN, 



de 



u = nx. 



(66) 



Relations express a well-known FS property of a solution: the product n'^\fl+T}?W, 
being one of invariants of the RG generator, is expressed in terms of some universal func- 
tion Mq of another invariant of the generator ()64j). and a form of this function is set by 
initial conditions, i.e. it is defined in terms of initial distribution functions of particles /q . 
Such representation has a common enough nature and is a typical property of BVP solu- 
tions obtained with the help of RGS which is pertinent to name ^-theorem (on analogy 
with known H-theorem). We briefly discuss this property in concluding remarks. 

In order to illustrate formulas (jU^ we apply them to plasma slab that contains ions 
of several, say two, types (the index g = 1, 2 corresponds to heavy and light ions, respec- 
tively) with initial Maxwellian velocity distribution functions, and the electrons obeying a 
two-temperature Maxwellian distribution function with densities and temperatures of the 
cold and hot components rico and Uho (^co + ^/io = Tlq Zqnqo) and Tc and Th, respectively. 
In this case the density distribution and, hence, the ion energy spectrum is expressed as 



Figure El demonstrates the plots of Afq for the following plasma parameters: Ti 2/Tc = 0.1; 
Th/Tc = 1000; nho/ZiUio = 5 x lO"'*. Block curves show dependence of a dimensionless 
"universal" density of plasma ions Nq = (nqo/nco)A/'q, referred to the maximal density 
of cold electrons, upon the dimensionless "coordinate" = {Js/LoY- Dotted curves 
present the distribution of the dimensionless density of cold and hot electrons (short and 
long strokes respectively). These curves demonstrate the high end of the energy spectrum 
for light ions {q = 2) with a sharp decrease, as well as the spatial separation of ion species 
which is of current interest in experiments on interaction of short laser pulses with thin 
foil targets (see, e.g., IS]). Similar results are obtained for more complex distribution 
functions jlH] and beyond the scope of the model used for the one-dimensional expansion, 
for example for spherically-symmetric expansion of a plasma bunch j45j . 

Summarizing the paragraph 14.31 let's note, that here a modified RG algorithm for 
calculation of RGS for non-local systems of the equations, and a procedure of prolongation 
of RG operators on solution functionals are used simultaneously. 




(67) 



where the function S is defined in the implicit form 




(68) 



9=1,2 



x{l + {ZqTnJmq))] - n/joexp [(1 - {Tc/Th))£] . 
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Figure 3: The curves describing the dependence of invariants of the RGS operator 
the "universal" density Nq of plasma ions - carbon ions (curves (C)) and protons (curves 
(H)) - is represented as a function of a dimensionless "coordinate" = LqY/ (l+fi^t^). 
For illustration here by dotted curves with short and long strokes the dependencies of a 
dimensionless density for hot and cold electrons upon are also shown. 

5 Conclusion 

The realization of the program that expands opportunities of application of the RG algo- 
rithm to problems of mathematical physics, is the goal of the paper, which is obviously 
formulated in section 1 for non-local problems. The appropriate class covers now (besides 
problems on a basis of DE, discussed in section 2 in an introductory example) the models 
containing non-local terms, including integral and inetgro-differential equations. 

This formulation preserves the former general scheme of construction of RG algorithm 
as four consecutive steps jl]. However the form of realization of these steps significantly 
varies, that is most brightly shown on first two stages of the algorithm (see figure Q), 
related to construction of non-local RG manifold and a calculation of admitted symmetry 
group. Here in view of absence of a regular computational algorithm (similar to Lie 
algorithm for DE) while performing of the second stage of algorithm various realizations 
are possible. As an illustration we choose and state in more detail the variant that is 
based on use of a canonical operator. 

In the following section 4 efficiency of use of procedure of prolongation the RG operator 
on non-local variables is shown with the purpose of the reduced description of the solution 
in terms of the integrated characteristic, the solution functional. It is essential, that the 
knowledge of a solution in an explicit form therewith is not required. 

The examples given here serve as an illustration of the program formulated in section 
3. They show application of new algorithm to integro-differential systems (expansion of a 
plasma bunch) and a procedure of prolongation on solution functionals (Cauchy problem 
for Hopf equation, evolution of a laser beam in nonlinear optics). The example presented 
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in 14. 31 is the most valuable, as here the calculation of RG symmetries for a solution of the 
non-local system of equations is supplemented by a procedure of prolongation of the RG 
operator obtained on the solution functional, the density of plasma particles, with the 
purpose of revealing its variation law. 

At the formulation and discussion of the analytical form of results an accent is made 
on a role of invariants of appropriate RG operators. The manifested common regularities 
are considered in the Appendix. 

The results of section 4 testify to universality of a method of RenormGroup Sym- 
metries. Therefore they allow to look forward to a further expansion of a class of the 
problems which can be investigated with the help of RGS method, and to new objects, 
for which the application of RG algorithm yet is not a standard procedure. 

Here we mean an infinite systems of looped integro-differential equations, similar to 
systems for correlation functions in statistical physics and to systems of the equations for 
generalized Green functions - propagators and vertex functions - in the quantum field 
theory. 

The work was partially supported by RFBR grant No. 05-0 1-00631, grant of Scientific 
School No.2339.2003.2 and ISTC project 2289. 
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6 Appendix. Invariant representation of solution and 
<l>-theorem 



In the main text of this paper the role of invariants of RG transformation in construction 
of BVP solutions was repeatedly marked. Here we consider the relation of representations 
of BVP solution with the concept of functional self-similarity and with the well-known 
principles of the group analysis. 

Let us remind that BVP solutions which are obtained with the use of RG algorithm 
are invariant solutions of RG operators. In group analysis of differential equations the 
explicit expression of solutions through invariants uses the well-known theorem of 
invariant representation of a regular (non-singular) manifold (see. [251 §18,ch.5], and 
|24j Vol.3, p.6]): 

Let a manifold M C admit a group Q. Suppose M be a nonsingular manifold 
of a group i.e., an infinitesimal operator of group G does not vanish identically 
on M . Then, M can be represented by a system of equations, left-hand sides of 
which are invariants of group Q, i.e., have the form: 

^k{Ji{z), J2{z), Jn-i{z)) = 0, k = l,...,s. (69) 

Here Ji(z), J2(-2), . . . , Jn-i{,z), z G M^, form a basis of invariant of Q. Hence, equations 
(j69j) with arbitrary functions oi N — 1 variables furnish the general form of non-singular 
invariant manifold of the group Q. In particular, it gives a transparent comment of the 
well-known H-theorem [lUj . 

Turn now to BVP solutions which are obtained with the use of RG algorithm. For RG 
invariant solutions there exists a more general statement as compared with H-theorem, 

^-theorem: An invariant solution of a boundary-value problem can be rep- 
resented by a system of equations of the form \69^) written down in terms of 
functional invariants (pi of the problem. 

These (pi are understood as invariants of appropriate functional transformations involving 
not only dependent and independent variables of the equations, but also parameters of 
boundary conditions, that is invariants of renormgroup transformations. 

In essence, $ -theorem is an analogue of the theorem of invariant representation with 
reference to solutions of BVP having the property of the functional self-similarity. In this 
case, one should consider sub-manifold TZA4 invariant with respect to renormgroup TZQ 
as a nonsingular manifold. 

In a special case, when in (|69|) s = 1 , and in the functional invariant (p{y, {a}) , 
containing the required function y , variables are separated, the solution can be written 
down in an explicit form close to the representation, which emerges from the H-theorem: 

y = </.-;($(...), {a}); </)=$(... ,0,,...), z = l,...,iV-l. (70) 

Here, the function is a reverse one to (p with respect to its first argument. Due to this 
solution y appears dependent not only on the remaining functional invariants, (pi , but 
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also on variables and parameters, {a} , entering into the invariant (f) . Thus, as well as for 
power self-similarity, BVP solutions are not, generally, invariants of RG transformations, 
but are expressed through certain combinations of invariants of RGS operators. 

Expressions (jUT|) for distribution functions of plasma particles in an expanding bunch 
serve as an example of BVP solution, being such invariants. 

As the second example, we take a QFT model with two coupling constants g and 
h. Here, invariant quantities, e.g., observed effective scattering cross-sections cr^{s) , are 
expressible in terms of RG invariants — two invariant coupling functions g{s/ fj,'^; g, h) , 
h{s/ fi^] g, h) and of the ratio rn? / s — by relations 

ay{s) =T.y[m^/s,g ,h) . (71) 

In turn, functions g and h should be found from system of two functional relations (see, 
e.g., eqs. (48.37) in Ref. [11]) 

G{y/x, g{x,y;g,h), h{x,y;g,h)) = G{y;g,h); 

H{y/x, g{x,y;g,h), h{x,y;g,h)) = H{y;g,h), (72) 

containing two arbitrary functions, G and H , of two arguments. Due to this, to find each 
of o"jy(s) , one needs to have explicit expressions for three defining functions T,,^, G and 
H. 

Note also that the procedure of numerical defining of the parameters g and h from 
boundary data (in fact, from observed quantities) involves at least two implicit relations 

(ED. 

At the same time, functionals of functions which determine, according to formulae 
(|65|) . the density distributions of particles of expanding plasma, are not invariants of 
the RG operator. Another example when a BVP solution is constructed with the help of 
invariants of RG transformations, but is not such invariant itself, is submitted by formulae 
fl50|) and for functionals in a problem of a beam refraction in a nonlinear medium. 
Generally, when it is impossible to express the BVP solution in an explicit form, one should 
use general formulae instead of the representation (fTO)) . 
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